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We investigate tiie kinetics of tlie A -\- B ^ Q reaction witli long-range attractive interaction 
V{r) ~ —r'^" between A and B or witli the drift velocity v ~ r'" in one dimension, where r is 
the closest distance between A and B. It is analytically show that the dynamical exponents for 
density of particles (p) and the size of domains [t) continuously vary with a when a < ac — /1/2, 
while that for the distance between adjacent opposite species {Iab) varies when a < cr^^ = 7/6. 
Beyond , diffusive motions dominate the kinetics, so that the dynamical behavior for diffusive 
systems is completely recovered. These anomalous behaviors with the two crossover values of a are 
supported by numerical simulations and the argument of effective repulsion between the opposite 
species domains. 
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The irreversible two-species reaction A + B ^ has 
been intensively and widely investigated as a basic model 
for various phenomena in physics 0, Q, chemistry 
and biology 0. When an A and a B particle meet on 
the same site, they instantly and irreversibly combine to 
form an inert species. Until now the studies for the re- 
action have been focused on understanding the effect of 
the fluctuations of the initial particle density or the global 
bias of particle motions on the kinetics without careful 
consideration of the interactions between A and B even 
for charged particles IllE0,llli,li,E3,Ii3- However, 
in situations where the kinetics of reaction is affected 
by certain attractive interactions between A and B such 
as the Coulomb interaction, the interactions should be 
much more important than the simple diffusive motions 
or global biases to understand the kinetics. Such sit- 
uations may include matter-antimatter annihilation in 
the universe, soliton-antisoliton recombination, charge 
recombination in clouds Ti] and electron-hole recombi- 
nation in irradiated semiconductor structures . 

It has recently been shown through a simple model 
that the underlying attractive interaction between oppo- 
site species leads to completely different scaling behav- 
iors from those studied so far In the model, the 
fluctuation-dominated kinetics leads to the segregation 
of alternating A— rich and _B— rich domains as usual f^. 
In addition, domain-boundary particles feel the attrac- 
tive interaction and are assumed to move to the adjacent 
opposite species domain with a constant drift velocity. 
Since the bias is a constant regardless of the distance be- 
tween A and B, the interaction strength in the model 
|l3l | is infinite. In contrast the particles inside the do- 
main are screened by the same neighboring particles and 
the motion of bulk particles is naturally assumed to be 
isotropic diffusion. As a result, the interaction causes the 
alternatively changing bias at domain boundaries which 
is neither the relative nor the uniform bias of Refs. 
and ^3 (Fig-df^))- In non-interacting systems, the den- 



sity decay has been known to depend on the motion and 
the mutual statistics of particles. For isotropic diffusions, 
the particle density scales as pU) ^ t^'^/^ in d di- 
mensions (d<4)|E|alllilS[l3- With the global 
relative drift, p{t) scales as p{t) ~ for d < 3 

0. With the uniform drift of both species, the hard- 
core (HC) constraint between identical particles leads to 
the scaling of p{t) ~ t~^/^ in one dimension With 
the infinite interaction, pit) scales as t~'^/^ regardless of 
the HC constraint jl,3l |. In one dimension, the uniform 
bias of HC particles and the infinite interaction lead to 
the same scaling law, but the scaling behaviors of ba- 
sic lengths are different. Hence the infinite interaction 
results in new dynamical scaling behaviors [l^ . 

The constant drift or the infinite attractive interaction 
is unnatural [l^ and cannot explain more general or real 
situations where the attractive interaction depends on 
the distance r between A and B. Some of physically re- 
alistic attractive interactions should be those described 
by a conserved attractive potential V{r) ~ — r"^*^. The 
drift velocity v is then given as v{r) ~ vTV^ ~ r~'^ . 
In this paper, we investigate the kinetics oi A + B ^ Q 
with the attractive potential V{r) ~ -^r^^"" or with the 
drift velocity to opposite species v{r) — r^" . In real- 
ity the motions of boundary particles cannot be deter- 
mined only by the interaction, because there should exist 
various noises, which are normally believed to make the 
background diffusive motions. In this reason, a bound- 
ary particle is assumed to stochastically move to its op- 
posite species domain with the rate p = (1 -I- i;(r))/2, 
and to its own domain with 1 — p for a unit time. Hence 
this model naturally includes the competition between 
the attractive interaction and the diffusive motion. The 
drifted motions at boundaries and the competition are 
expected to lead rich and interesting scaling behavior as 
a varies. Fig.n(b) shows the space-time trajectories for 
a = 0.3. Using the trajectories schematically depicted 
in Fig. n] (c), we analytically derive asymptotic scaling 
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FIG. 1: (Color online) Snapshot of trajectories of A + B 
with the attractive interaction between A and B of (a) cr = 
and (b) a — 0.3. (c) The magnified schematic trajectories 
of adjacent opposite species domains. Subscripts {1,2,. ..,n} 
indicate the order of the positions of particles from a given 
domain boundary. 



behavior. Intriguingly, we find that all dynamical ex- 
ponents for various densities and lengths continuously 
vary when a < ac- Furthermore the scaling behaviors 
for the present model manifest unique anomalous behav- 
iors, which have not seen in dynamical critical models 
with the continuous varying exponents |l4|. In addition, 
the anomalous behavior is that there exist two different 
crossover values of ac- cFc for the density of particles (p) 
and the size of domains {£) is CTc ~ 1/2, while that for 
the distance between adjacent opposite species (Iab) is 
a^^ = 7/6. Hence the interaction completely changes 
the scaling behavior of kinetics when a < crf^. How- 
ever, when (7 > cr^^, the kinetics for diffusive motions is 
completely recovered 0. As in the a = case, the HC 
constraint is irrelevant due to the isotropic diffusions in- 
side domains. We also numerically confirm our analytical 
results. 

With the equal initial density pa(0) — Pb{0), parti- 
cles are randomly distributed on an one-dimensional lat- 
tice of size L. When the selected particle has two op- 
posite species neighbors such as {A - ■ ■ A - ■ ■ B), the cen- 
tral particle (A) hops to the opposite species {B) with 
rate p = {1 + v{r))/2 and to the same species (A) with 
rate 1 — p, where v{r) = r^'^ with a > and r is the 
distance between A and B. Otherwise, particles diffuse 
isotropically. In the region of a length £, the number of A 
species is initially Na = /Oa(0)£± \/ pa{0)£ and the same 
for Nb- After a time t ^ particles travel throughout 
the whole of the region, and annihilate by pairs. The 
residual particle number is the number fluctuation in the 
region, so we have the relation Na ~ or ~ l/v^ 
for a given length ^ |H ■ As the processes evolve, the 
system becomes a collection of alternating A-rich and B- 
rich domains. To characterize the structure of segregated 
domains, we introduce three length scales as in Ref. 1^. 



The length of the domain {£) is defined as the distance 
between the first particles of adjacent opposite species 
domains |^. The length Iab is defined as the distance 
between two adjacent particles of opposite species, while 
£aa{£bb) is the distance between adjacent A{B) particles 
in a A{B) domain. These lengths scale asymptotically as 
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A bulk particle inside single species domains diffuses 
isotropically and the mean position of the bulk particle 
is not substantially changed until it becomes a bound- 
ary particle. Once it becomes boundary, it drifts to 
its neighboring opposite species with velocity vlr) un- 
til it annihilates in the midway between two opposite 
species domains. From this situation in mind we now 
analytically calculate the scaling behaviors of the kinet- 
ics when the drift velocity is v{r) ^ r~'^ . For an AB 
pair with the distance £ab, the time interval for anni- 
hilation of the AB pair is tab — £ab/v{£ab) ^ ^"ab ■ 
For cr > 0, the space-time trajectory of a particle is 
arch-shaped. In the tr = case the trajectory is very 
close to a pentagon 0. These arch-shaped trajectories 
should be self-similar (self-afRne) fractal structures, be- 
cause they should have the scaling symmetry due to the 
power-law scaling behavior (1). A typical base unit of 
the self-similar arch-shaped trajectories of adjacent op- 
posite domains are schematically depicted in Fig. ^c). 
This base unit allows us to calculate a time ti needed to 
remove the unit of size £ surrounded by one scale larger 
ones. Then the size of the larger unit increases by £ dur- 
ing Ti and we have 



M/dt r^£/Tl, 



(2) 



which gives the dynamic exponent z. In following cal- 
culations, we consider only the mean positions of bulk 
particles, and assume £AA{t) to be a constant during the 
annihilation of the base unit. After a smaller unit is 
completely annihilated, the remainder of particles redis- 
tribute over the larger unit increased by the size of the 
annihilated unit. Hence we approximate £AAit) = ■ ■ ■ — 
tAA{t + Tn) = • • • = £AA{t + Tp) during the annihilation 
of a smaller unit. 

As only boundary particles of each domain have two 
opposite species neighbors, the annihilation of boundary 
particles comes from the attractive interaction. It takes 
a time ti = £'^j^g for the first boundary pair, Ai and 
Bi in Fig. ^ (c), to annihilate. The second pair {A2 
and B2) from the boundary isotropically diffuse during 
time Ti until Ai and Bi annihilate. After the time ri, 
the second pair becomes a new boundary pair, and the 
^2 ~ {£ab + ^aaY^^"^ is needed for the annihilation. So 
it takes time T2 ~ ti -I- t2 in total for the second pair 
to annihilate. Similarly, the nth pair from the initial 
boundary will annihilate after r„ ~ r„_i -t-i„, where tn = 
{£ab + {n — 1)£aa)'^'^^ for n > 2. From the recurrence 
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relation of r„, we find 

(3) 

- ^A^B + [i^AB + nlAAY+^ ~ ^'abM^AA ■ 

The second line is obtained by integrating out the sum- 
mation over k. As the number of particles Ni and the 
length Iaa in ^ domain of size t scale as Nf^ ^ ^/i and 
£aa ~ 1/p Vl, the time T£ needed to annihilate the 
domain of size £ in the base unit is given by 

Ti - iAl' + mB + er+' - rA^y^^^ ■ (4) 

Because of ^ > £ab by definition, we finally find the 
leading scaling of ti as 

r,^r+3/2 , (5) 

and z is given as z = a + 3/2 from Eq. ||2Jl. Intriguingly, 
z increases continuously and linearly with a. However 
z cannot increase beyond the upper bound Zc — 2, be- 
cause domains cannot spread more slowly than those of 
random diffusion. Hence for a > (Tc = 1/2, we have 
z = 2 regardless of a. The scaling of £aa is easily ob- 
tained from the relation £aa ~ V£ ~ t^^^-*^. One finds 
zaa = 2z = 2a + 3. Since £aa with the diffusion only 
scales as £aa ^ t^^'^ the critical value of a is also 
af^ = 1/2. For £ab, we consider the change of density 
during time tab |S{ • During tab , one pair of AB particle 
annihilate only between boundaries, the change of par- 
ticle density is given as dp/dt ~ —Pab/tab, where pab 
is the density of AB pairs.. Using relations, p ~ 1/v^, 
PAB = 1/^, and TAB - £aV, one finds £ab - t^^"^^^ 
with zab — z = a + 3/2. Interestingly, £ab follows the 
scaling of £. However, £ab with the diffusion only scales 
as £ab ^ t^'^ 9] and thus a^. of £ab is cr^^ ==7/6. The 
scaling exponents of three lengths and their critical val- 
ues of tr, beyond which the diffusive scaling behaviors 
recover, are summarized as follows. 

z = a + 3/2 , ZAA = 2z , zab = z , , . 

(Tc = 1/2 , a^^ = 1/2 , = 7/6 . ^""^ 

From the scaling of lengths, asymptotic decays of var- 
ious densities can be extracted. The density of to- 
tal particles (p = 2/?^), adjacent pairs of same species 
{paa = Pbb) and adjacent pairs of opposite species 
{pab) scale as 

P^t-", PAA - t""^^, pab ~ t""^^ . (7) 

ks p ^ l/\f£, we have p ~ i^^/^^ with a = l/2z = 
1/(2(7 + 3). PAA follows the same scaling of p due to 
PAA ~ 1/£aa ~ l/Vi so PAA t~°' with aAA = a. 
Finally pab scales as pab ~ 1/^, which leads to pab ~ 
t~i/^ with aAB — 1/z ~ !/("' + 3/2). As densities with 
the diffusion only scale as p ~ paa ^ t^^^^ and pab ~ 



the upper bound of cr is CTc = 1/2 for all densities. 
All decay exponents of various densities are simply given 
as 

a = l/2z, aAA^a, uab^^/z , (8) 

and the critical value of ct is Cc = 1 /2 for all densities. For 
cr = 0, the all scaling behaviors of the constant drift are 
fully recovered (l3|. As shown in Eqs. (6) and (7), there 
exists two different crossover values of cr, Cc (= cr^^ — 1/2 
and a^^{— 7/6), where Cc < ^c^- It is very pecu- 
liar that there exists two crossover values. However the 
inequality ctj, < directly comes from the anoma- 

lous scaling of £ab with diffusion only. With diffusion 
only, adjacent opposite domains are effectively repulsive 
by preferential annihilations of nearby AB pairs. Hence, 
opposite species pairs are further apart than the typical 
interparticle distance_r4A and, as a result, £ab increases 
anomalously as t^^^ |9| . This anomalous scaling of £ab 
leads to another crossover in addition to the crossover by 
diffusion. 

As the strength of the interaction becomes weak, the 
effect of the random fluctuation by diffusion becomes 
strong. The competition between the fluctuation by dif- 
fusion and the drift by the interaction leads to the con- 
tinuously decaying exponents and the critical values of 
a above which diffusive motions dominate the kinetics. 
On the other hand, since the interaction maintains the 
Galilean invariance of the domain structure, the predic- 
tions © and (|SJ) are expected to be independent of the 
HC constraint. 

To confirm our analytic results numerically, we now 
present the simulation results. With pa{0) = Pb{0), A 
and B particles distribute randomly on a chain of size 
L. In the simulations we consider both HC particles and 
the particles without the HC constraint, which we call 
bosonic particles. In the model with HC particles there 
can be at most one particle of a given species on a site. 
In the bosonic model there can be many identical parti- 
cles on a site. All the simulations are done on the chain 
of size L = 3 X 10^ and with pa{0) = Pb{0) = 0.1. 
We average densities and lengths over 100 independent 
runs. Fig. 12 (a) shows densities (inset) and their effective 
exponents defined as —a{t) — log[p(2i)/p(t)]/log2 for 
a — 0.1 case of HC particles. We estimate a = 0.317(5), 
ctAA — 0.3125(25), and uab = 0.625(25) respectively. 
All results agree well with the prediction (jSJ ; a = uaa = 
0.3125 and uab = 0.625 for a = 0.1. Fig. El (b) shows 
the various lengths and their effective exponents defined 
asl/z(t) = \og[£{2t)/£{t)]/log2 for a = 0.1. We estimate 
z = 1.60(5), ZAA = 3.20(5) and zab = 1.60(5) which 
also agree very well with the prediction ©; z = 1.6, 
Zaa — 3.2 and zab = 1-6 for a — 0.1. For bosonic parti- 
cles, we also obtain nearly the same results. 

To check the crossover from the continuously varying 
dynamic scaling to the completely diffusive one, we per- 
form Monte Carlo simulations under the same initial con- 
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FIG. 2: (Color online) Effective exponents of densities (a) 
and lengths (b) of the model with HC constraint for a = 0.1. 
In each panel, two horizontal lines from top to bottom show 
the predictions (a) a = 0.313 and a^s = 0.63 (b) 1/ z = 0.625 
and 1/ zaa = 0.313. 
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FIG. 3: (Color online) Effective exponents of densities (a) 
and lengths (b) of the model with HC constraint for a — 1.25. 
In each panel, horizontal lines from top to bottom show the 
prediction (a) a = oaa ~ 1/4 and uab ~ 1/2 (b) 1/z = 1/2, 
1/zAB = 3/8 and 1/zaa = 1/4. 



ditions for a = 1.25(> a^^ > CTc)- Fig. shows simu- 
lation results of HC particles for a — 1.25. We estimate 
a = 0.25(1), aAA = 0.24(1) and a^s = 0.47(1), which 
agree reasonably well with the values without the drift in 
one dimension; a = aAA = 1/4 and aAB = 1/2 9]. For 
lengths, we estimate 1/z — 0.50(5), 1/zaa = 0.26(1), 
and 1/zAB = 0.375(5) which also agree well with the 
values without the drift, 1/z = 1/2, 1/zaa = 1/4 and 

l/zAB^5/8 0. 

We plot the estimates of all exponents versus a and 
the lines for the predictions © and Q both for HC par- 
ticles and for bosonic particles in Fig. 4. All exponents 
except 1/zAB continuously vary along the predicted lines 
for (T < CTc = 1/2. In contrast 1/zab varies along the pre- 
dicted hne for a < a^^ = 7/6. Beyond CTc (or <J^^), each 
exponent takes the value of the diffusive system with- 
out the interaction regardless of a. The exponents for 
HC particles are nearly identical to those for the bosonic 
particles, and this result comes from the irrelevance of 
the HC constraint to the kinetics due to the isotropic 
diffusion inside domains. 

In conclusion, we investigate the kinetics of irreversible 
reaction A + B ^ with the attractive interaction 




FIG. 4: (Color online) Plot of dynamical exponents versus 
a. Main (inset) plot show the exponents of HC (bosonic) 
particles. Solid lines in each panel from top to bottom 
correspond to the predictions; 1/zab = 1/ (a -I- 3/2) and 
1/zaa = 1/(2(7-1-3). Horizontal dashed lines from top to bot- 
tom correspond to l/z = 1/2, 1/zab ~ 3/8, and 1/zaa = 1/4 
respectively. 



V{r) ^ — j-^2cr |-,g^-^yggjj ^ Qj- ■y^rit}l the drift veloc- 

ity v{r) ^ r^'^ to opposite species at segregated domain 
boundaries. The drift leads to arch-shaped space-time 
trajectories of particles with which we analytically de- 
rive asymptotic scaling behaviors, and numerically con- 
firm them. Intriguingly, the drift results in continuously 
varying scaling behavior by certain critical values, Uc or 
. above which diffusive motions dominate the kinet- 
ics. In our derivations, we only consider the mean posi- 
tions of particles and neglect the expansion of domains by 
the random fluctuations during the annihilations of do- 
mains. Hence when the interaction dominate the kinet- 
ics, the fluctuation does not affect the asymptotic scaling 
behavior. 
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